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"t^ ■ Abstract 

^^ , We consider the problem of distributed formation control of a large number of vehicles. An individual vehicle in the formation 

^^ is assumed to be a fully actuated point mass. A distributed control law is examined: the control action on an individual vehicle 

depends on (i) its own velocity and (ii) the relative position measurements with a small subset of vehicles (neighbors) in the 
^^ formation. The neighbors are defined according to an information graph. 

.^^ , In this paper we describe a methodology for modeling, analysis, and distributed control design of such vehicular formations 

whose information graph is a D-dimensional lattice. The modeling relies on an approximation based on a partial differential 
equation (PDE) that describes the spatio-temporal evolution of position errors in the formation. The analysis and control design 
is based on the PDE model. We deduce asymptotic formulae for the closed-loop stability margin (absolute value of the real 
part of the least stable eigenvalue) of the controlled formation. The stability margin is shown to approach as the number of 
^ ■ vehicles A'^ — >■ oo. The exponent on the scaling law for the stability margin is influenced by the dimension and the structure of 

the information graph. We show that the scaling law can be improved by employing a higher dimensional information graph. 

in 

CO ' Apart from analysis, the PDE model is used for a mistuning-based design of control gains to maximize the stability margin. 

Mistuning here refers to small perturbation of control gains from their nominal symmetric values. We show that the mistuned 
design can have a significantly better stability margin even with a small amount of perturbation. The results of the analysis with 



o, 

(^ , the PDE model are corroborated with numerical computation of eigenvalues with the state-space model of the formation. 



I. Introduction 

We consider the problem of controlling a group of vehicles so that they maintain a desired formation geometry while 
following a desired trajectory. The desired formation geometry is specified in terms of desired relative positions between pairs 
of vehicles. The desired trajectory of the formation is supplied to a subset of the vehicles, which are called reference vehicles. 
The problem is relevant to a number of applications such as formation control of aerial, ground, and autonomous vehicles 



IM- 



In many of these applications, a centralized control 



for transportation, surveillance, reconnaissance and mine-sweeping 

solution that requires all-to-all or all-to-one communication is impractical. This motivates distributed control architectures 

where an individual vehicle exchanges information only with a small set of other vehicles to make control decisions. 

Each vehicle is modeled as a fully actuated point mass. This means that (i) the dynamics of each coordinate of the vehicle's 
position are modeled using a double integrator, (ii) the coordinate dynamics are decoupled, and (iii) an independent force 
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actuation is assumed for each coordinate. A distributed control law is examined: the control action on an individual vehicle 
depends on (i) its own velocity and (ii) the relative position measurements with a small subset of vehicles (neighbors) in the 
formation. The neighbor relationship is defined according to an information graph, which has been recognized to play an 
important role in closed-loop stability of the formation ||5i|6[]. A node in the graph represents a vehicle, and an edge connecting 
two nodes represents direct information exchange between those nodes. This information exchange may occur due to one 
vehicle measuring the other vehicle's relative position by on-board sensors, or due to one vehicle communicating its state 
information to the other via a communication channel. The information graph is undirected, meaning that if vehicle i can get 
information about vehicle j, then j can get information about i. 

The objective of this paper is to study how the stability margin (the absolute value of the real part of the least stable 
eigenvalue) of the closed-loop scales with the number of vehicles, structure of the D-dimensional information graph, and the 
choice of the control gains. For a specific case, when D = I, the stability margin of the platoons and its dependence of 
asymmetry in control gains was examined in our own earlier work [7]. The extension to 2D formations appears in [8J. This 
paper is an extension of these previous works. 

In this paper, we restrict ourselves to information graphs that belong to the class of _D-dimensional (finite) lattices. A formal 
definition of lattice appears in Section [III see Figure |2] for a few examples. Lattices arise naturally as information graphs when 
the vehicles in the group are arranged in a regular pattern in space and the exchange of information occurs between pairs 
of vehicles that are physically close. However, lattices also allow for a flexibility to model much more general information 
exchange architectures. They are often used as information graph in concensus and vehicular formation problems 191]. In this 
paper, we make an important distinction between the dimension of the position vector of a vehicle and the dimension of the 
lattice that defines the information graph. For example, a one-dimensional platoon may have a two-dimensional lattice as its 
information graph (see Figure [U. 

A. Related literature 

An analysis of the stability margin is important to understand the scalability of control solutions as the number of vehicles 
in the formation, N, increases. In the formation control literature, the scalability question has been investigated primarily for 
a one-dimensional vehicle formation, which is usually referred to as a platoon. An extensive literature exists on the platoon 



control problem; see 11101 - 11411 and references therein. The most widely studied information exchange structures for distributed 
control of platoons are predecessor following control and bidirectional control. In predecessor following control, every vehicle 
uses information from the vehicle immediately ahead. In bidirectional control, each vehicle uses information from the vehicle 
immediately ahead of it and the one behind it. Scenarios in which information exchange occurs with vehicles beyond those 
physically closest, are studied in 115111611 . The focus of much of the research in this area has been on the so-called symmetric 
control, in which every vehicle uses the same control law. Such a simplifying assumption is motivated in part by a lack of 
tools for analysis and design of distributed control laws. The symmetry assumption is used to simplify the design and analysis. 



References that studied non-symmetric control design include 116 . 



IJ. 



For platoons, the distributed control architectures with symmetric control are known to scale poorly, both in terms of 
closed-loop stability margin and sensitivity to external disturbances. In a symmetric bidirectional architecture, the least stable 
closed-loop eigenvalue approaches zero as A^ increases (TD. This progressive loss of stability margin causes the closed-loop 
performance to become arbitrarily sluggish as the number of vehicles, N, increases. Small stability margin can also lead to long 



transients due to initial conditions, which can resuh in control saturation III8II . It is worthwhile to point out that the stability 
margin for a platoon is known to scale poorly as a function of N even with the centralized LQR control il9ll . In addition 
to the loss of stability margin, the sensitivity of the closed-loop platoon to external disturbances increases without bound 



2111 or slinky-type effect H13I1 . String instability is 



2311 ■ Non-symmetric control design 



as a function of N. This effect is also referred to as string instability 020 , 

observed for both symmetric predecessor following and symmetric bidirectional control H22i 

within the bidirectional architecture was proposed in JtI] that helps improve the closed-loop stability margin. A non-symmetric 

I I 

control design within the framework of predecessor architecture was proposed in [17], which ameliorates string instability at 
the expense of control gains that increase without bound as N increases. 

Control of platoons with inter-vehicle communication that allows for information exchange with vehicles that are not just 
nearest neighbors was considered in llisl II6II . It was concluded in 11511 that to eUminate string instability with symmetric 
control, the number of vehicles that each vehicle communicates with has to grow without bound as N increases. It was shown 
in il6ll that heterogeneity in control gains does not significantly alter string instability if certain constraints are imposed on 
integral absolute error and high frequency response of the loop transfer function. 

Bamieh et. al. studied controlled vehicle formations with a _D-dimensional torus as the information graph (|9|]. Scaling 
laws with symmetric control are obtained for certain performance measures that quantify the sensitivity of the closed-loop to 



stochastic disturbance. It is shown in 



m 



that the scaling of these performance measures with N is strongly dependent on 



the dimension D of the information graph. In 112411 . Pant et. al. introduced the notion of mesh-stability for two-dimensional 
formations with a "look-ahead" information exchange structure, which refers to a particular kind of directed information flow. 
The scenario considered in our paper, with undirected information graphs, does not fall under the look-ahead information 
exchange structure. 

B. Contributions of this paper 

In this paper we describe a methodology for modeling, analysis, and distributed control design of vehicular formations 
whose information graph belongs to the class of /^-dimensional lattices. The approach is to use a partial differential equation 
(PDE) based continuous approximation of the (spatially) discrete platoon dynamics. Just as a PDE can be discretized using a 
finite difference approximation, we can carry out the procedure in reverse: the spatial difference terms in the discrete model 
are approximated by spatial derivatives. The resulting PDE yields the original set of ordinary differential equations upon 
discretization. This approach is motivated by earlier work on PDE modeling of one-dimensional platoons DTI]. The PDE model 
is used for analysis of stability margin and for mistuning-based design of distributed control laws. 

There are two contributions of this work that are summarized below. 

First, we obtain scaling laws of the stability margin of the closed-loop formation with symmetric control. We show that 

the stability margin scales as 0{-^) where ni is the number of vehicles along a certain axis of the information graph. By 

"1 

choosing the structure of the information graph in such a way that rii increases slowly in relation to N, the reduction of the 
stability margin as a function of N can be slowed down. In fact, by holding ni to be a constant independent of the number 
of vehicles N, the stability margin can be bounded away from zero even as the number of vehicles increase without bound. 
It turns out, however, that keeping rii fixed while N increases causes the number of reference vehicles to increase. When the 
information graph is a square _D-dimensional lattice (equal number of nodes on each side of the lattice), the stability margin 
scales as 0{ ^2^0 ) in the limit of large TV. This formula is a generalization of the estimate given in |7] for a one-dimensional 
formation. 



The second contribution of this work is a procedure to design asymmetric control gains so that the stability margin scaling 
law is significantly improved over that with symmetric control. For the case of square information graphs, we show that an 
arbitrarily small asymmetry in the proportional control gains from their nominal symmetric values results in stability margin 
scaling as 0{ „}/o )■ In contrast to the 0{jAjd) scaUng seen in the symmetric case, this is an order of magnitude improvement. 
The resulting control design is called a mistuning-hcLsed design since the control gains are changed only slightly, i.e., mistuned, 
from their values in the nominal, symmetric case. Mistuning-based approaches have been used for stability augmentation in 
several applications 125142811 . and recently for distributed control of one-dimensional platoons ||7|]. 

The advantage of using a PDE-based analysis is that the PDE reveals, better than the state-space model does, the mechanism 
of loss of stability and suggests the mistuning-based approach to ameliorate it. Numerical computations of eigenvalues of the 
state-space model of the formation is used to confirm the scaling laws with symmetric as well as mistuned control. Although 
the PDE model approximates the (spatially) discrete formation dynamics in the limit N — > oo, numerical calculations show 
that the conclusions drawn from the PDE-based analysis holds even for small number of vehicles. 

The remainder of this paper is organized as follows. Section |ll] presents the problem statement and the main results of this 
paper Section|III]describes the state-space and PDE models of the formation control problem. Analysis and control design results 
together with their numerical verification appear in Sections |IV] and [V] respectively. In Section |Vl] we present time-domain 
simulations to illustrate these results, and comment on various aspects of the proposed design and analysis methodology. 

II. Problem statement and main results 

A. Problem statement 

We consider the formation control of N identical vehicles. The position of each vehicle is a Z?s -dimensional vector (with 
-Ds = 1, 2 or 3); Ds is referred to as the spatial dimension of the formation. Let p^ £ M be the d-th coordinate of the i-th 
vehicle's position, whose dynamics are modeled by a double integrator: 

pf^=uf\ d=l,...,i?„ (1) 

where u\ G M is the control input (acceleration or deceleration command). The underlying assumption is that each of the 
Ds coordinates of a vehicle's position can be independently actuated. We say that the vehicles are fully actuated. The spatial 
dimension Dg is 1 for a platoon of vehicles moving in a straight line, Dg — 2 for a formation of ground vehicles and Dg = 3 
for a formation of aerial vehicles flying in the three dimensional space. 

The control objective is to make the group of vehicles track a pre-specified desired trajectory while maintaining a desired 
formation geometry. The desired formation geometry is specified by a desired relative position vector A^ j :— p*{t) ~p*{t) for 
every pair of vehicles (i, j), where p*{t) is the desired trajectory of the vehicle i. The desired inter-vehicular spacings have to 
be specified in a mutually consistent fashion, i.e. A^ j = A^^ + A^.j for every triple i, j, k. Desired trajectory of the formation 
is specified in the form of a few fictitious "reference vehicles", each of which perfectly tracks its own desired trajectory. The 
reference vehicles are generalization of the fictitious leader and follower vehicles in one-dimensional platoons [TI llOl ll9|]. A 
subset of vehicles can measure their relative positions with respect to the reference vehicles, and these measurements are used 
in computing their control actions. In this way, desired trajectory information of the formation is specified only to a subset of 
the vehicles in the group. In this paper we consider the desired trajectory of the formation to be of a constant-velocity type, 
so that Ajj's don't change with time. 

Next we define an information graph that makes it convenient to describe distributed control architectures. 
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(a) The desired formation geometry of a ID spatial platoon (b) The desired formation geometry of a 2D (c) The information graph for both the ID 
with 6 vehicles and 3 reference vehicles. spatial vehicle formation with 6 vehicles and platoon and the 2D formation shown in (a) 

3 reference vehicles. and (b). 

Fig. 1. (a, b): Two distinct spatial formations that have the same associated information graph (c). Red (filled) circles represent reference vehicles and black 
(unfilled) circles represent real vehicles. Dashed lines (in (a), (b)) represent desired relative positions, while solid lines represent edges in the information 
graph. 



Definition 1: An information graph is an undirected graph G — (V, E), where the set of nodes V = {1, 2, . . . , A^, A^ + 
1, . . . ,iV + Nr} consists of A^ real vehicles and N^ reference vehicles. The set of edges E C V x V specify which pairs 
of nodes (vehicles) are allowed to exchange information to compute their local control actions. Two nodes i and j are called 
neighbors if {i,j) G E, and the set of neighbors of i are denoted by A/i. D 

Note that information exchange may or may not involve an explicit communication network. For example, if vehicle i 
measures the relative position of vehicle j with respect to itself by using a radar and uses that information to compute its 
control action, we consider it as "information exchange" between i and j. If a vehicle i has access to desired trajectory 
information then there is an edge between i and a reference vehicle. 

In this paper we consider the following distributed control law, whereby the control action at a vehicle depends on i) its 
own velocity and ii) the relative position measurements with its neighbors in the information graph: 



M) - 



E -^..^-^f" - ^) - ^^ " -*''')^ ^ = 1, . . . , AT, 



(2) 



jeM 



where «*('*) is the d-th component of the desired velocity of the formation, k,J is the proportional gain and 6; s is the 
derivative gain. Note that all the variables in (|2]i are scalars. It is assumed that vehicle i knows its own neighbors (the set Mi), 
desired spacing A^- \ and the desired velocity v*'^'^\ 

Example 1: Consider the two formations shown in Figure [T] (a) and (b). Their spatial dimensions are Ds — 1 and Dg — 2, 
respectively. The information graph, however, is the same in both cases: 

V = {1,2,...,9}, E = {(1,2), (1,4), (1,7), (2, 3), (2, 5), (2, 8), (3, 6), (3, 9), (4, 5), (5, 6), (7, 8), (8, 9)}. 

A drawing of the information graph appears in Figure [T](c). Although the information graph is the same, the desired spacings 

Ai j's are different in the two formations. For example, Aj 5 7^ in the one-dimensional formation shown in Figure [T] (a) 

whereas A2 5 = in the two-dimensional formation shown in Figure [T](b). 

In this paper we restrict ourselves to a specific class of information graph, namely a finite rectangular lattice: 

Definition 2 (D-dimensional lattice): A _D-dimensional lattice, specifically a ni x n2 x • • • x njj lattice, is a graph with 

nin2 ■ ■ ■ TiEi nodes. In the D-dimensional space MP , the coordinate of i-th node is i :— [ii, . . . , id]^ , where ii E {0, 1, . . . , {ni — 
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(a) A ID 4 lattice. 
Fig. 2. Examples of ID, 2D and 3D lattices. 
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(c) A 3D 2 X 3 X 3 lattice. 



1)}, 12 € {0, 1, • • ■ , (?^2 — 1)}, • • • and ijj G {0, 1, ... , {ud — 1)}. An edge exists between two nodes i and j if and only if 
\\i — j\\ = 1, where || • || is the Euclidean norm in M^. A rti x rt2 x • • • x n^ lattice is denoted by Z„jxn2x-xnr>- With a 
sUght abuse of notation, "the i-th node" is used to denote the node on the lattice with coordinate i. D 

Figure |2] depicts three examples of lattices. A /^-dimensional lattice is drawn in M^ with a Cartesian reference frame whose 
axes are denoted by xi,X2, ■ ■ ■ , xd- Note that these coordinate axes may not be related to the coordinate axes in the physical 
space R^^ 

In this paper an information graph G is always a lattice Z„jxri2-xni5, where nin2 ■ ■ .riD = N + Nr- For a given N, 
the choice of Nr,D,ni,n2, ■ ■ ■ ,nD serves to determine the specific choice of the information graph within the class. An 
information graph is said to be square if rii — 1 = ?t.2 = . . . = no- 

For the ease of exposition and notational simplicity, we make the following two assumptions regarding the reference vehicles 
and the distributed control architecture (|2|: 



^(d) 



(d) 



Assumption 1: For each (i, j) G E, the gain k,^ .^ does not depend on d, and for each i G V, h\ does not depend on d 



D 



Assumption 2: The reference vehicles are arranged so that a node i in the information graph corresponds to a reference 
vehicle if and only ifzi=ni — 1. D 

Assumption [T] means that the local control gains do not explicitly depend upon the coordinate d. Such an assumption is not 
restrictive because of the fully actuated assumption. If the local control gains are allowed to depend upon d then one could 
repeat the analysis of this paper separately for each value of d. Note that the assumption does not mean that the control gains 
are spatially homogeneous; for example, the control gains fcL ^ 7^ ^a ■) for 'he same {i,j) G E. 

Assumption |2] means that all reference vehicles are assumed to be arranged on a single "face" of the lattice, and every 
vehicle on this face is a reference vehicle. Assumption |2] implies that N = (ni — 1)712 • ■ -no and Nr = 712 ■ • "£>■ Other 
arrangements of reference vehicles do not significantly change the main conclusions of this paper. Some of these extensions 
are discussed in Sec. |Vl] 

As a result of the Assumption [l] we can rewrite Q as 



Ui = ^ -ki^,^j){p, -pj - A,,j) - bi{p, - V*), 



(3) 



where the superscript (d) has been suppressed. 



Remark 1: The dimension D of the information graph is distinct from the spatial dimension Ds- Figure [T] shows an example 
of two formations in space, one with Dg = 1 and the other with Z?, = 2. The information graph for both the formations is 
the same 3x3 two-dimensional lattice, i.e., D — 2. On account of the fully actuated dynamics and Assumption [T] the spatial 
dimension Dg plays no role in the results of this paper The dimension of the information graph D, on the other hand, will 
be shown to play a crucial role. 

Remark 2: Analysis of the control law (O is relevant even when there are additional dynamic elements in the controller 
There are several reasons for this. First, a dynamic controller cannot have a zero at the origin, for that will result in a pole- 



zero cancellation, causing the steady-state errors to grow without bound as N increases ll23ll . Second, a dynamic controller 



cannot have an integrator either If it does, the closed-loop platoon dynamics become unstable for a sufficiently large value 
of A^ i23ll . Thus, any allowable dynamic compensator must essentially act as a static gain at low frequencies. The results 
of ilsL I23II indicate that the low frequency behavior is the dominant factor in the control of large networks of agents with 
double integrator dynamics. Hence, the issues that arise with the control law ^ are also relevant to the case where additional 
dynamic elements appear in the control law. 

B. Main result I: Stability margin with symmetric control and D-dimensional information graph 

Definition 3: The stability margin is the absolute value of the real part of the least stable eigenvalue of the closed-loop 
system. D 

Definition 4: The control law (O is symmetric if all the vehicles use the same control gains: fc(i _,) = ko, for all {i,j) E E 
and bi = 5o for all i e V, where fco and bo are positive constants. 

The first main result gives an asymptotic formula for controlled formation with symmetric control: 

Theorem 1: Consider an A^-vehicle formation with vehicle dynamics ^ and control law (|2|l, with Assumptions [T] and |2l 
With symmetric control, the stability margin of the closed-loop is given by the formula 

that holds when ni -^ cxi. D 

We remark that the stability margin depends only upon rii - the number of vehicles along the xi axis. The xi axis is special 
because it is normal to the face with the reference vehicles; see Assumption|2] In the PDE model, the boundary condition is of 
the Dirichlet type on this face (see (l28Tl). Analogous estimates also hold with different arrangement of the reference vehicles 
(see Section [Vl] for details). 

a) Square information graph: For a square information graph, TV = {ni — l)?i2 . . . n£, — {ni — 1)^, and we have the 
following corollary: 

Corollary 1: Consider an A^-vehicle formation with vehicle dynamics ([T]l and control law (|2]i, with Assumptions [T] and |2] 
When the information graph is a square D-dimensional lattice, the closed-loop stability margin with symmetric control is given 
by the asymptotic formula 



li 



The special case of Corollary [T] for D = 1 was established in 

The result from Corollary [T] shows that for a constant choice of symmetric control gains fcp ^nd b^, the stability margin 
approaches as A^ ^ 00. The dimension D of the information graph determines the scaling. Specifically, the stabiUty margin 



scales as 0(l/iV^) for ID information graph, as 0{1/N) for 2D information graph, and as 0{1/ N'^/'^) for 3D information 
graph. Thus, for the same control gains, increasing the dimension of the information graph improves the stability margin 
significantly. In practice, this may require a communication network with long range connections in the physical space. Note 
that an information graph is only a drawing of the connectivity. A neighbor in the information graph need not be physically 
close. 

Remark 3: It was shown in [19] that the closed-loop stability margin for a circular platoon approaches zero as 0{l/N^) even 
with the centralized LQR controller. It is interesting to note that distributed control (with an information graph of dimension 
D > 1) yields a better scaling law for the stability margin than centralized LQR control. 

b) Non-square information graph: It follows from Theorem [T] that by choosing the structure of the information graph 
in such a way that ni increases slowly in relation to N, the loss of the stability margin as a function of N can be slowed 
down. In fact, when ni is held at a constant value independent of N, it follows from Theorem [T] that the stability margin 
is a constant independent of the total number of vehicles. More generally, consider an information graph with rii = 0{N'^), 
where c e [0, 1] is a fixed constant. Using Theorem [T] it follows that S — 0{1/N'^^) as A^ — > cx). If c < -^, the resulting 
reduction of S with N is slower than that obtained for a square lattice; cf. Corollary [T] This shows that within the class of 
D dimensional lattices (for a fixed D), certain information graphs provide better scaling of the stability margin than others. 
The price one pays for improving stability margin by reducing rii is an increase in the number of reference vehicles. This is 
because the number of reference vehicles iV^ is related to ni by Nr = N/{ni — 1) (see Assumption |2l). 

It is important to stress that not all non-square graphs are advantageous. For example, if ni ~ 0{N) and n2 through no 
are 0(1), it follows from Theorem [T] that the stability margin is S* = 0(1/A^^). This is the same trend as in a 1-D information 
graph. In this case, we can say that the D dimensional information graph effectively behaves as a one dimensional graph. 

Figure [3] shows a few examples of information graph that are relevant to the discussion above. The 2D information graph 
shown in Figure [5] (a) has ni = 0(1) and 712 = 0{N), whereas the one in Figure [3] (b) has ni = 0{N) and n2 = 0(1). The 
graph shown in Figure [3] (c) is approximately square, both rii and ri2 are 0{^/N). 

Figure |4] provides numerical corroboration of the discussion above. The stability margin as a function of N for three distinct 
2D information graphs (that are described in Figure [3]) are shown in this figure. The stability margin is computed by computing 
the eigenvalues of the closed-loop state matrix; the state space model is described in (fT2b in Section |III1 The control gains 
used are ko — 0.01, bo — 0.5. The plots show that the formula Q in Theorem [T| makes an excellent prediction of the trend 
of stability margin. The asymptotic nature of the result in Theorem [T| (and Corollary [TJ is seen from the plot: the prediction 
becomes more and more accurate as A^ increases. 

C. Main result II: Stability margin with non-symmetric control and D-dimensional information graph 

The second main result of this work is that for a fixed information graph, the scaling law for stability margin can be improved 
by choosing a non-symmetric control law. We call the resulting design a mistuning-hased design because it relies on small 
changes from the symmetric control. The improvement is achieved by making small perturbations to the proportional gains 
alone k(ijy, changing derivative gains alone do not have the same disproportionate effect, it only effects the 0{l/n1) term in 
the stability margin. The mistuning-based design and the resulting scaling law is summarized with the aid of the following 
theorem: 
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(a) Non-square information grapli, S = 0(1) 



(b) Non-square information graph, S = 0(1/Af^) (c) "Approximately" square information graph, 

S = 0{l/N) 



Fig. 3. (a) A 2-dimensional information graph in which the first dimension is held constant, resulting in a stability margin that is independent of N, 
S = 0(1). (b) A 2-dimensional information graph that is "asymptotically" ID (as N —>■ oo) since the size of the first dimension increases linearly with TV, 
resulting in a stability margin scaling law S = 0(1/N'^), which is the same as that with an ID information graph, (c) A 2-dimensional information graph 
in which both sides are of length 0{y/N), for which we have S = 0(1/N), the same behavior as that of a square 2D graph. 



Theorem 2: Consider an A^-vehicle formation with vehicle dynamics ([T]) and control law (|2]i under Assumptions [T] and |2] 
with nominal symmetric control gains /cq and fep- Now consider the problem of maximizing the stability margin by designing the 
proportional control gains fc(,j j), where the gains are required to satisfy |fc(i j) — fcol < e for every {i,i) G E, with e E (0, fco) 
being an arbitrary and small pre-specified constant. For vanishingly small values of e, the optimal control gains of the i-th 
vehicle (i — 1, . . . , N) are given by: 



«(i,ii 



fco + £, fc(i.ii-) — ko — s, k(^ij) = fco for all other neighbors j, 



(6) 



where i^^ denotes i's neighbor in the positive xi direction (in the drawing of the information graph) relative to node i and 
i^~ denotes i's neighbor in the negative xi direction. The resulting stability margin is given by 
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bo ni — 1 n 

The formula is asymptotic in the sense that it holds when tii, . . . , no -^ oo and e — > 

For the special case of a square information graph, we have the following corollary. 

Corollary 2: For a vehicular formation of N vehicles with square information graph and mistuned control design described 
in Theorem |2] the stability margin is given by 
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l/D + ^^ at2/D>' 



where e is defined in Theorem |2] 



(8) 

D 



We note that the additional information needed by each vehicle i to implement the mistuned control comprises of (i) the 
parameter e and (ii) the knowledge of which one of its neighbors is the neighbor i^+ and which is i^^. The special case of 
Corollary |2] for D = I was established in ||7|]. 

Comparing Theorems [T] and |2] (similarly, Corollaries [T] and |2]l, we see that the effect of mistuning is to introduce a square 
root in the stability margin formula. Thus, even for a small e, mistuning can improve the closed-loop stability margin by a 
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Fig. 4. Stability mai'gin predicted by Theorem [T] for a vehicle formation with information graphs of vaiious "shapes" as shown in Figure [1] The legend 
"SSM" means computed from the "state space model" U2t . which is presented in Section ITlTl For the first case, ni — 1 = 5 and n2 = N/5. Theorem [T] 
predicts that in this case S = 0(1) even as N ^ oo. In the second case, 712 = 5 and ni — 1 = N/5, which leads to 5 = 0{1/N'^). The third case is that 
of a square information graph, ni — 1 = 712 = vJV, which leads to 5 = 0(1/N). Theorem [T] and corollary [T] predicts the stabihty margin quite accurately 
in each of the cases. The control gains used in all the calculations are fco = 0.01 and bo = 0.5. 

large amount, especially when N is large. Numerical verification of the conclusion of Theorem |2] is presented in Section [Vj 
see, in particular. Figure [TO] Figure |5] depicts the optimal mistuned control gains for the case where the information graph is 
a 3 X 3 lattice. 



III. Closed-loop dynamics; State-space and PDE models 

A. State-space model of the controlled vehicle formation 

The dynamics of the i-th vehicle is obtained by combining the open loop dynamics ([T]i with the control law ([3]l, which 
yields 



Pi = X! ~hhj)iPi - Pj - Aij) - bi{pi - V*), 1 = 1,, 



,N. 



(9) 



Let pI {t) denote the desired trajectory of the i-th vehicle. The trajectory is uniquely determined from the trajectories of the 
reference vehicles and the desired formation geometry. For example, suppose the trajectory of a reference vehicle r is v*t. If 
the d-th coordinate of the desired gap between a vehicle i and the reference vehicle r is A^ .^, then the d-th coordinate of the 
desired trajectory of i is p*'^'^\t) = v*^'^H + A] ,^. 

To facilitate analysis, we define the following coordinate transformation: 



p,; := Pt - _p, 



Pi =Pi-V 



Substituting (fTOl i into (|9]l, we have 



Pt^ ^ -fc(j,j)fe -Pj) - hp^- 



(10) 



(11) 
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Fig. 5. Optimal mistimed proportional control gains for a formation of 6 vehicles with 3 reference vehicles whose information graph is a 3 X 3 lattice. In 
general, each vehicle in the 2-dimensional information graph has 4 proportional gains (as shown on the top of the right figure), k,^ ii+)> ^(i ii-). ^^(i i2+) 



and k 



(i,«^-) 



The gain k 



(i,il + ) 



is the proportional gain with respect to the neighbor in the positive xi direction of vehicle i. The other three proportional 



gains have similar interpretations. Notice that the optimal mistuned control gains are achieved by simply making k,^ ^i+-. larger than the nominal gain fco by 
e and making k,^ ^i-% smaller than the nominal gain kg by e. The other proportional gains remain the same as the nominal gain. 



Since the trajectory of a reference vehicle is assumed to be equal to its desired trajectory, pi = if i is a reference vehicle. 
To express the closed-loop dynamics of the formation compactly, we define: 

P := [■P1,P2, ■ ■ ■ ,PNf, V := p = [■pi,P2, ■ ■ ■ ,PNf 
Using (fTTl i. the state-space model of the vehicle formation can now be written compactly as; 



^ ^ = Alp 



(12) 



where ip := [p; v] is the state vector and A the closed-loop state matrix. 

Example 2 (Example 1 contd.): Consider the ID and the 2D spatial formations depicted in Figure[T|(a) and (b), respectively. 
The information graph for both these formations is the same and drawn in Figure [TJc). We will now show that the closed-loop 
dynamics of both the formations are the same; cf. Remark [T] Specifically, let us examine the dynamics ^ for the vehicle 
i — 2. For the ID formation (Ds — 1), we have 



p'^' 



u(i) f„(i) 

"(2,1)1-^2 



rf' 



.(1)^ 



.(1) ^Jl) 



■P^^ 



(1)^ 



^2,l7 '^(2,3)^/^2 ^"3 ^2,3 



'^'(2,5) yi'2 



Pf' 



A^!l) 



(13) 



For the purpose of illustration, we focus on the third term on the right hand side of the above equation, and note that the 
desired trajectories are defined with respect to reference vehicle 7 (it can be defined with respect to any reference vehicle): 

*(1) _„*(l)|j. a(i) 



Using the notation in Eq. ( fTOl l. the third term in the right hand side of ( fT3] l can now be expressed as 



(14) 
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(1) lA^) 

(2,5)v/'2 



-Pi^^ 



.(1) 



.(1) ^^(1) 



*(1) 



.(1) 



-kyj.,{p^^>-p^^^), 



"(2,5) ^ 



^2,5 J — '^(2,b)^P2 +P2 Pb P5 ^2,5^ 

where the first equality follows from ( fT4b and A2 7 — Ag 7 — A2 5, which follows from the definition A^ ._, = p* 
evaluating the other terms in a similar manner, we obtain 
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Fig. 6. A pictorial representation of the i-th vehicle and its four nearby neighbors in a 2D information graph. i^+ stands for the neighbor of the i-th vehicle 
in the xi positive direction relative to vehicle i, and i^^ stands for the neighbor of the i-th vehicle in the xi negative direction relative to vehicle i. And 
4-^+ and i^~ can be interpreted in the same way. 



In case of the formation with spatial dimension Dg — 2, we examine the dynamics of the second component of the position 
vector of vehicle 2: 

-kZ,{p'i'-p'^'-^'il)^b'i\i^^' --*''')■ (16) 

For this formation, the desired trajectories are also defined with respect to reference vehicle 7, 

p*W^„*W^ + ^W^ p*W^„*(fc)^ + ^W^ fc^l^2, (17) 

so that the third term on the right hand side of (fTST i can be expressed as 

,(2) , (2) (2) .(2)^ ,(2) ,~{2) , *(2) ~(2) *(2) .(2)^ ,(2) , ~{2) ~(2)^ 

^2,5)(P2 -Pi ~ Ay) = -k\^[z){P2 +P2 -P5 ~P5 ~ ^U) = -fc(2.5)(P2 " ^5 )' 

where the second equality follows from ([TtI i and A2 y — Ag y == Aj 5, which follows from the definition Aij = p* — p*. 
Repeating this procedure for each of the terms, one obtains: 

"(2) ,(2) .-(2) ~(2)n ,(2) ,~{2) ~(2)x ,(2) /~(2) ~(2)x ,(2) / ~(2) ~(2)n , (2) ~(2) ,, p. 

P2 = -fc(2,l)(P2 -Pi )- fc(2,3)(P2 -P3 ) - fc(2,5)(P2 " Ps ) " ^(2,8)(P2 " Ps ) " ^2 ^2 • (18) 

Under Assumption [T| that the gains are independent of d, ( fTsT i has the same structure as (fTsT l. The same holds for all the 
vehicles, which shows that the closed-loop dynamics (fT2] | depends only on the information graph. D 

Our goal is to analyze the closed-loop stability margin with increasing number of vehicles N and to devise ways to improve 
it by appropriately choosing the controller gains. While in principle this can be done by numerically computing the eigenvalues 
of the matrix A, such a computation does not clearly reveal the dependence of stability margin on N, control gains, graph 
structure etc. For this purpose, we approximate the dynamics of the spatially discrete formation by a partial differential equation 
(PDE) model that is valid for large values of N. The PDF model is used for analysis and control design. 

B. PDE model of the controlled vehicle formation 

For a given choice of the information graph, the i-th vehicle has the coordinate i — [ii, J2, . . . , in]^ in R^- We interpret pi 
as a function of the coordinate i. In the following, we consider a continuous approximation of this function to write a PDE 
model. 
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Fig. 7. Original lattice, its redrawn lattice and a continuous approximation of a discrete value function defined on the redrawn lattice, (a) is a 2D information 
graph for a formation with 3x3 vehicles and 3 reference vehicles, (b) shows a redrawn information graph of (a), so that it lies in the unit 2-cell [0, 1]'^. 
(c) gives a pictorial representation of continuous approximation of a discrete function whose values are well defined on the nodes in the redrawn lattice as 
shown in (b). 



For the i-th node with coordinate i — [ii, . . . , io]'^, we use 1"^+ and i'^~ to denote the nodes with coordinates [ii, . . . , id-i, id+ 
1, id+i, ■ ■ ■ , io]'^ and [ii, . . . , id-i,id — 1, id+i, ■ ■ ■ , *d]"^, respectively. For D = 2, a node i in the interior of the graph and 
its four neighbors, i.e., i^+, i^~,i^+, and i^^, are shown in Figure |6] The dynamics ( fTTT i can now be expressed as: 

D D 

Pi = -^ k{i^id+){pi - p,d+) - ^ k^,id--^{pi - p,d-) - bip,, (19) 



We define, 



kf'f^' ■.=k 



(1,^+) 



k 



{i.i''-) 



d=l 



Kj ' : — fc(jjd+') K^jjd--), 



de{l,...,D}. 



where the superscripts / and b denote fmnt and back, respectively. Substituting (|20l ) into ( fT9] ). we have 

-{Pi -Pi^-) 



Pi 



D j^dj+b^^dj^b 

hp^ = - ^ ^ ^-^ (p. - p,a+) - Y^ 

d=l d=l 



D_ j^dj+b _ ^dj-b 



(20) 



(21) 



2 '^' "' ^ ^ 2 

To proceed further, we first redraw the information graph in such a way so that it always lies in the unit _D-cell [0, 1]^, 



irrespective of the number of vehicles. Note that in graph-theoretic terms, a graph is defined only in terms of its node and edge 



sets. A drawing of a graph in an Euclidean space, also called an embedding 12911 . is merely a convenient visualization tool 



For the rest of this section, we will consider the following drawing (embedding) of the lattice Z„j 



x---xrio 



in the Euclidean 



space M . The EucUdean coordinate of the i-th node, whose "original" Euclidean position was [ii, . . . , id] , is now drawn at 
position [iici, i2C2, . . . , idCdV , where 



1 



Cd ■ = - 



d = l,. 



,D. 



(22) 



nd-l 

Figure |7] shows an example, where the original lattice, shown in Figure |7] (a), is redrawn to fit into [0, 1]^, which is shown in 
Figure |7](b). 

The starting point for the PDF derivation is to consider a function p{x,t) : [0, 1]^ x [0, oo) — > M defined over the unit 
D-cell in R^ that satisfies: 



Pi(t) ^ p{x,t)\s=li,ci 



J2C2,...,»r)C_Dj 



(23) 
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Figure He) pietorially depiets the approach: funetions that are defined at discrete points (the vertices of the lattice drawn in 
[0, 1]^) will be approximated by functions that are defined everywhere in [0, 1]^. The original functions are thought of as 
samples of their continuous approximations. We formally introduce the following scalar functions k^, k^, b : [0, 1]^ -^ M (for 
d e {!,..., D}) defined according to the stipulation: 



,12C2,---,IDCd\ ' 



bi — b{x)\g^[i-^ci,i2C2,...AnCoV ■ 



(24) 



In addition, we define functions k^ , n,^ 



^+' k^-'' : [0, 1]^ 



ki+\x):^kf,{x) + k'M, 
Due to JTM . these satisfy 

To obtain a PDE model from (l2Tl i. we first rewrite it as 

D 

Pi + biPi =y^fc; 



ki-\x):^ki{x)-k',{x). 



(25) 



k 



dj-b _ U-b 



k^ [X)\x=liici 



l2C2....,loCn\ 



d=l 



d.f-b {Pi'i+ - Pid-) 
Cd ^ 



E 

d=l 



kf'-^'^'' 2 iPi''+ - 2p» + Pt'i- ) 



(26) 



and then use the following finite difference approximations for every d E {1, . . . , D}: 



Pid+ - Pid- 



2cd 
p^d+ — 2pi + Pid- 



dp{x, t) 



dxd 
d^p{x, t) 



dxd^ 



X = lliCi,l2C2,...,loCo\ 



X = [liCi,l2C2,...,lDCn\ 



We emphasize that xi,. . . ,xd above are the coordinate directions in the Euclidean space in which the information graph is 
drawn, which are unrelated to the coordinate axes of the EucUdean space that the vehicles physically occupy. Substituting the 
expression (|22] | for Cd, ( l26] l is seen as a finite difference approximation of the following PDE: 



V ut ut / 



Kx)jAp{x,t) = Y, 



p{x,t), 



(27) 



^j^ V n-d - 1 dxd 2{nd - 1) dxd^ 
The boundary conditions of PDE jZTT i depend on the arrangement of reference vehicles in the information graph. If there are 
reference vehicles on the boundary, the boundary condition is of Dirichlet type. If there are no reference vehicles, the boundary 
condition is of the Neumann type. Under Assumption |2] the boundary conditions are of the Dirichlet type on that face of the 
unit cell where the reference vehicles are, and Neumann on all other faces: 



dp 

p{l,X2,...,XD,t) == 0, {0,X2,...,XD,t) = 0, 

oxi 



dp 
dxd 



(28) 



(f,i)=0, f = [xi,...,Xd_i,0 or l,a;rf+i,...,XD]^, (d > 1). 



If other arrangements of reference vehicles are used, the boundary conditions may be different. For future use, we rewrite the 
PDE dag as 

where the linear operator C is defined as 



^—^b{x)i-)pix,t)^£{--, 



dt 



dxd ' dxd^ 



)p{x,t), 



(29) 



d d 



^I^T.'' fl-r ,2 •' ■ 2^ 



D U-b 



ki-\x) d , fcr (^) d 



f+b. 



' dxd ' dxd 



nd~ I dxd 2{nd - 1)^ dxd'^ 



(30) 



It can be verified in a straightforward manner that the PDE dZTJ i yields the original set of coupled ODEs iT% upon discretization. 
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IV. Stability margin with symmetric control 
A. PDE-based analysis of stability margin 

Recall that in case of symmetric control we have 

k(ij)=ko, V(«,j)eE, bi^ba, Vi G V, 

where fco and 60 are positive scalars. In this case, using the notation in ( |20] l and (l24l l. we have 

k^/\x) = 2ko, fc^-''(f)=0, bix) = bo, d=l,...,D. 



The PDE (|29l l simplifies to a damped wave equation: 



(£+^4)^(^^'*)=^°(aS^)^(^^'^)' ^''^ 



where £o{x) is the Laplacian operator: 



where 



g2 q2 q2 q2 

o-l - (^^^ 1 ly ' d^l,...,D, (33) 

are the wave-speeds. The closed-loop eigenvalues of the PDE model require consideration of the boundary value problem 

'^o(a^)'^(^) = -^^i^)^ (34) 

For the given boundary condition of ( l28T l, the eigenvalues (different from the eigenvalue of PDE) and eigenfunctions of £9 
are respectively given by 

'{III ~ Ij'^'X 2 I /; \2 2 I t f-i \2 2 



2, / (2/i - 1)' ?i ?j) \ 

- "" '^oUCm - 1)2 + K - 1)2 + ■ ■ ■ + (n,, - 1)2;- 



4>f{x) — cos ( ) cos(/27ra;2) • • • cos{Ibttxd)- (35) 

where /i G {1, 2, . . . } and I2, ■ ■ ■ ,Id G {0, 1,2,...}. We use the notation / = (/i, . . . , Id) to denote the wave vector and Ap 
(j>i{x) to denote the associated eigenvalue and eigenfunction given by (|35] |. After taking a Laplace transform of both sides of 
the PDE ( I3TI 1 with respect to t, we get Vrj{x, s) = where T-" := s^ + &os — >Co and rj{x, s) —^ ip^{x)ai{s) is the Laplace 
transform of p{x,t) with aj^s) being its weights. Note that (f)^- is also the l-th basis of the null space of operator V. The 
eigenvalues of the PDE turn out to be the roots of the characteristic equation: 

s^ + bos + X^=0, (36) 

where s as the Laplace variable and A,- is an eigenvalue of (|34] i. The two roots of (|36T l are 



-bo ± Jbl - A\^ 
sf := \ . (37) 

We call si the l-th pair of eigenvalues. If the discriminant in ( l37b is positive, both the eigenvalues are real-valued. In this 
case, si is closer to the origin than sZ; so we call si the ^th less-stable eigenvalue. The least stable eigenvalue is the one 
among them that is closest to the imaginary axis, and the stabihty margin is the absolute value of its real part: 

Smin = minst, S := |-Re(smin)|- (38) 

r ' 
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Provided each of the n^'s are large so that the PDE ( |27] i with the boundary condition (|28] l is an accurate approximation of 
the (spatially) discrete formation dynamics (fT2] i under Assumption |2] the least stable eigenvalue of the PDE (|29] l provides 
information on the stability margin (see Definition |3]l of the closed-loop formation dynamics. We are now ready to prove the 
Theorem [U that was stated in Section HIl 

Proof of Theorem [7] Consider the eigenvalue problem for PDE ( [3T| i with mixed Dirichlet and Neumann boundary condi- 
tions ( |28] |. Since the less stable eigenvalues are given by st = |(— 60 + \/^o ^ ^-^f)- ^^ ^^ discriminant 6§ — 4A,-is positive. 



both of the eigenvalues are real-valued. In this case, st is closer to the origin than sZ\ so we call st the Z-th less-stable 
eigenvalue. It follows from ( |37| | that the least stable among them is the one that is obtained by minimizing A^-over the Z?-tuples 
(Zi, . . . , lo). Using dSSl l. this minimum is achieved at Zi = 1, Z2 = • • • = /d = 0, where A(l, 0, . . . , 0) = 0.257r^fco/(ni — 1)^. 
Therefore, 

.2; ) 



Smin = mm S^ = -r-(-l + 1 - TT? ■-, x^ 



TT^fcn 



46o(ril - 1)2 






where the last equality holds when ni ^ 1- 
from the equation above. 



T^yko 



Due to the definition of stability margin ( I38I I. the result follows immediately 
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Fig. 8. Numerical comparison of closed-loop eigenvalues with symmetric control predicted by the state-space model (SSM) J12t and PDE model (39) with 
mixed Dirichlet-Neumann boundary conditions. Eigenvalues shown are for a square information graph with 26 X 25 nodes (625 vehicles), and the control 
gains used are fco = 0.01, bo = 0.05. Only a few eigenvalues are compared in the figure. PDE eigenvalues are computed using a Galerkin method 13011 . 



B. Numerical comparison of eigenvalues between SSM and PDE 

We now present numerical computations that corroborates the PDE-based analysis. We consider a 26 x 25 square two- 
dimensional information graph with symmetric control. The gains are 

fc(i,ii+) = fc(j,ii-) = fc(i,i2+) = k(^i,i2-) = ko ^ 0.01, 6j = 6o = 0.05. 
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The associated PDE model is given by 

92 



(^ + 6„|Mf,t).( 



fco 



92 



fco 



92 



(ni-l)2 9a;2 {n^ - lY dxi 



p{x,t). 



(39) 



The eigenvalues of the state matrix A in ( fT2] l are compared against the eigenvalues of the PDE ( |39] l with mixed Neumann- 
Dirichlet boundary conditions in Figure [8] The eigenvalues of the PDE are computed numerically using a Galerkin method 
with Fourier basis 13011 . The comparison in Figure |8] shows that the PDE eigenvalues match the state-space model eigenvalues 
well, especially the ones close to the imaginary axis. Figure |9] shows, as a function of N, the stability margin computed from 
the PDE and the state-space model. The prediction from the asymptotic formula (|5]l in Corollary [T] is also shown. We see from 
Figure |9] that the least stable eigenvalue of the closed-loop is well captured by both the PDE model as well as the asymptotic 
formula © that is derived from analysis of the PDE. 
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Fig. 9. The stability margin of the closed-loop formation dynamics with symmetric control (fcg = 0.01 and bo = 0.5) as a function of number of vehicles: 
the legends of SSM, PDE and Corollary [T] stand for the stability margin computed from the state space model, from the PDE model, and from the asymptotic 
formula {5) in Corollary [T] 



V. MiSTUNING-BASED CONTROL DESIGN 

With symmetric control, one obtains an 0(-^) scaling law for the least stable eigenvalue because the coefficient of the 

"1 

g|-^ term in the PDE dSTT l is 0{-^) and the coefficient of the g|- term is 0. Any asymmetry between the forward and the 
backward gains will lead to non-zero fc;J^ (x) and the presence of 0(;j-) term as coefficient of g|-. By a judicious choice 
of asymmetry, there is thus a potential to improve the stability margin from O(^) to 0( — ). The subsequent analysis shows 
that this is indeed so, and a control design is proposed to achieve the 0( — ) trend. One should also note that this insight into 
the control design problem is difficult to obtain from the examination of the state matrix A. 



A. Reducing loss of stability by mistuning 

In this section, we consider the problem of designing the control gain functions kj^{x) and k^{x) so as to improve the 
stability margin over symmetric control. Specifically, we consider the problem of minimizing the least-stable eigenvalue Smin 
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of the PDE ( |29] ) by changing the control gains slightly (mistuned) from their values in the symmetric case. We begin by 
considering the forward and backward position feedback gain profiles 

k{^{x) ^ko + eU^ix), k'^^ix) ^ko + efc^(x), (40) 

where e > is a small parameter signifying the amount of mistuning and k^{x), fc^(af) are functions defined over [0, 1]^ that 
capture gain perturbation from the nominal value fco- Define 

fc2(f) := yjx) + fc^(f), ~kT{x) := ~k!,{x) - fc^(f). (41) 

Due to the definition of kj^ and kj^ in dZSl l. we have 

kf/\x) = 2fco + efc^(x), k^d'^'i^) = £C(^)- 

The mistuned version of the PDE (|29T l is thus given by 

We study the problem of improving the stability margin by judicious choice of k^ (x) and fc™ (x) while keeping the gains k^ (x) 
and fc^(x) within certain pre-specified bounds. The results of our investigation, described in the following sections, provide a 
systematic framework for designing control gains in the formation by introducing small changes to the symmetric design. 

To design the "mistuning" profiles k^{x) and fc™(a?') to minimize the least stable eigenvalue Smin, we first obtain an explicit 
asymptotic formula for the eigenvalues when e is small. The result is presented in the following theorem. The proof appears 
in the Appendix. 

Theorem 3: Consider the eigenvalue problem of the mistuned PDE (|42] | with mixed Dirichlet and Neumann boundary 
condition (l28T l. The least stable eigenvalue is given by the following formula that is valid when e — !• and rii, n2, ■ ■ ■ , njj -^ oo: 

smin = 4°i„ - ^T^TT^ — TT / ^™(^) ^i" (^^i) '^^i - ^JTT^ — T^ / ^i(^) cos2( Jxi) dx^ + 0{e^), (43) 

where s^^l^ is the least stable eigenvalue without mistuning, i.e., of PDE (|3TT i with the same boundary conditions. D 

It follows from Theorem |3] that to minimize the least stable eigenvalue, one needs to choose only fc™(x) carefully; all other 
fc™'s and all fc|'s can be set to 0. The reason is that only k™{x) and fcf(a;) affect the least stable eigenvalue, and the term 
involving fcj(x) is of order l/(rii — 1)^, whereas the term involving k'^{x) is of order l/(ni — 1). For large ni the effect of 
the function fc™(a?') on the least stable eigenvalue will be far greater than that of fcf (£"). Therefore, we choose 

k'a{x) = = k^{x) for d = 2,...,D, and fc^(f) = 0. 

This means that the perturbations to the "front" and "back" gains satisfy k^ (x) — k\{x) = for d = 2, . . . ,D. For d — 1, the 
choice kl{x) = leads to 

k({x) = -klix) ^ k'^ix) = 2k({x). 

The most beneficial gains can now be readily obtained from Theorem|3] To minimize the least stable eigenvalue with kf{x) = 0, 
we should choose fc™(a^) to make the integral J^ fc™(x) sin(7rxi)(ia;i as large as possible, which is achieved by setting k™{x) 
to be the largest possible value everywhere in the unit cell. This result is summarized in the next Corollary. 
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Corollary 3: Consider the problem of minimizing the least-stable eigenvalue of the PDE ( |421 ) with mixed Dirichlet and 
Neumann boundary condition ( |28] ) in the limit as e — ?> by choosing ki{x), k\{x) e L°°([0, 1]) with the constraint that 
11^1 (2;)||oo — \\ki{x)\\ao = 1, where || • ||oo denotes the sup-norm. The solution to this optimization problem is given by 

fcf (f) = 1, fcj(f) = -1, Vxe[0,l]^- □ 

The proof of Theorem |2] now follows in a straightforward manner from Corollary |3] 

Proof of Theorem^ Note that ensuring |fc(ij) — fco I < e in the formation is equivalent to keeping |fc;^—fco| <£and|fc^— fcol <e 
for d = 1, . . . , _D in the PDE domain; cf. ( |24] |. This is equivalent to keeping ||fc;J ||oo < 1 and ||fc^||oo < 1 for each d\ cf. ( |40] |. 
In this case, the optimal gains are those given in Corollary [5] It follows from (l24l l that the optimal gains for the vehicles are 

k{^,^l+) = (fco + efcf (^))U=[»ici,...,»dcd]t = fco + £, V i G V 

k(i,ii-) = (fco + ek\{x))\3=[i^ci.,...,iDCDV = fco - £, V z G V 

fc(i.i''+) — fcfi.i''-) — ko, (i>l,VzGV. 

The resulting least stable eigenvalue is, from Theorem [5] 



/•l 9 1 

— — / sin(7ra;i)dxi + sfj^^ = -e— — + 0( — ), 

bo[ni - 1) Jo oo(ni - 1) nf 

since s^^jj^ — 0{l/nl). The result follows upon taking absolute value of Smin- 



'%. 



B. Comparison of eigenvalues between mistuned SSM and PDE 

Figure [TO] depicts the numerically obtained mistuned and nominal eigenvalues for both the PDE and state-space model for 
a 2D square information graph. The nominal control gains are fco — 0.01, 6o — 0.5, and the mistuned gains used are the ones 
shown in Figure |5] with e = 0.001. The figure shows that 

1) the closed-loop poles match the PDE eigenvalues accurately over a range of N; 

2) the mistuned eigenvalues show large improvement over the nominal case even though the controller gains differ from 
their nominal values only by ±10%. The improvement is particularly noticeable for large values of N, while being 
significant even for small values of N . 

For comparison, the figure also depicts the asymptotic eigenvalue formula given in Theorem |2] The improvement in the 
stabiUty margin with mistuning is remarkable since the gains are changed from their symmetric values by only ±10%. Another 
interesting aspect of the result in Corollary |2] is that the improvement from 0{1/N^^^) to 0{1/N^^^) can be achieved by 
arbitrarily small changes to the nominal gains. In addition, the optimal mistuned gain profile is quite simple to implement. For 
a vehicle formation with arbitrary dimensional information graph and with a maximum variation of ±10% from the symmetric 
gains, the optimal gains are obtained by letting fc(i ^1+) be 10 percent larger than the nominal gain fco and letting fc(i,ii-) be 
10 percent smaller than the nominal gain. 

VI. Discussion 

A. Relationship between the stability margins of the coupled-ODE and PDE models 

In this paper, all the analysis and control design are based on the stability margin of the PDE model, which is an approximation 
of the coupled-ODE model under the assumption that each ni (i G {1,2,..., _D}) is very large. This raises the question: how 
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Fig. 10. Stability margin improvement by mistuning for a vehicle formation with 2D square information graph. The nominal control gains are fco = 0.01, 
bo = 0.5, and the mistuned gains used are the ones shown in Figure [s] with e = 0.001. The symbol N on the x-axis is the number of vehicles and S on 
the j;-axis is the stability margin. The legends "Nominal SSM" and "Nominal PDE" stand for the stability margin computed from the state-space model and 
the PDE model, respectively, with symmetric control. The legends "Mistuned SSM" and "Mistuned PDE" stand for the stability margin computed from the 
state-space model and PDE model, respectively, with mistuned control. We see that the (i) the PDE model predicts the stabiUty margin quite accurately, and 
(ii) the StabiUty margin is improved significantly by mistuning control design even with ±10% variation from the symmetric gains, especially for large N. 



large is the difference between the stability margin of the PDE (continuous problem) and the coupled-ODE (discrete problem) 
model? In this section, we provide an analysis on the difference between the stability margins of the continuous and the discrete 
problems, which we call the stability margin approximation error. The results are summarized in the following lemma. 

Lemma 1: Consider an A^-vehicle formation with vehicle dynamics ([1]) and control law (|2]i, under Assumptions [1] and |2] 
With symmetric control (respectively, mistuning design), the stability margin approximation error between the PDE model ( |27l ) 
with boundary condition ( l28T l and the discrete model is 0(l/nf ) (respectively, 0{l/n1) + 0{e^)). D 

In particular, for a square information graph, the stability margin approximation error bounds for symmetric control and 
mistuning design are 0(l/iV^/^) and 0{l/N'^^^) + 0{e'^) respectively. 

Recall that for symmetric control (respectively, mistuning design), the stability margin scales as 0(l/n^) (respectively, 
0(1/711)). Comparing with the above lemma, we can see that the PDE model provides an accurate approximation to the 
coupled-ODE model, and the approximation error can be ignored even for a moderate value of 711, which is the number of 
vehicles along the xi axis of the information graph. For the ease of description, we only provide the proof for a formation 
with 1-dimensional information graph, i.e. the case D = I. Figure [TT] depicts a picture of the ID information graph. The proof 
for higher dimensional case follows in a similar manner, upon using the closed form expressions of the eigenvalues for the 
discrete case 03 ill . 

Proof of Lemma U] Given D ^ I, consider the following coupled-ODE and PDE models: 

Pi + biPt = -k(iA-i){Pi -Pi-i) - fc(i,i+i)(pj ~Pi+i), 
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Fig. 11. Pictorial representation of a ID information graph. 

d'^p{xi,t) dp{xi,t) ^ k(~''{xi)dp{xi,t) fcf+''(xi) d'^p{xi,t) 

9t2 "*■ ^^'^ ai " m - 1 a^i 2(ni - 1)' dxi^ ■ 

Using the optimal control gains given in ^ of Theorem ID the above models are simplified to: 

Pi + baPi = -(fco - e)(Pi - Pi-i) - (fco + £)iPi - Pi+i), (44) 

d'^p{xi,t) dp{xi,t) ^ 2 dp{xi,t) fco 02p(xi,t) 

9t2 +''0 at "^m-l 9a;i ^(n, -1)2 fei2 ■ ^ ^^ 

Notice that when e = 0, it corresponds to the symmetric control case. Now consider the following discrete and continuous 

eigenvalue problem: 



1) 



^Xpi = -(fco -e)(Pi -Pi-i) - (fco +£)iPi ~Pi+i), 



(46) 



where i £ {1, 2, . . . , rii — 2}, and for the 0-th vehicle, there is no neighbor behind it, so its equation is given by —Xpo 
— (fco + s){po — pi). And also, recall that the reference vehicle indexed by "ni — 1" has the property that Pm-i — 0. 

2 dp{x,t) fco d'^p{x,t) 



2) 



-fj,p{x,t) — e 



ni-1 dxi (m - 1)^ dxi"^ 



where the boundary condition is given by g^' ' = 0, p(l,t) = 0. 

For the discrete eigenvalue problem, we can write it compactly as Ap = Lp, where p := [po,P2, • 

as follows: 

fc'O + S —fco — £ 

-fco + e 2fco 



,Pni-2 



(47) 



F and L is defined 



-fco 



-fco 



(48) 



[32]. For the case of 



2fco —fco — £ 

-fco -I- £ 2fco 
For the symmetric control case {e = 0), the least eigenvalue of matrix L is given by 4fco sin^ 2(2n -i) 
mistuning design, under the assumption that e is small, we can use matrix perturbation method to compute the least eigenvalue 
of L (see [33]). Combining the results, we have the least eigenvalue for the discrete eigenvalue problem: 

2(1 -fees 2;^) 



A = 4fcosin „_ XI- 

2(2ni - 1) 2ni - 1 

By Taylor series expansion theorem, the above eigenvalue can be expressed as 

fco7r2 fco7r2 2 1 



+ 0(£^). 



A 



0{e ) + higher order terms. 



4(ni - 1)2 4(ni - 1)3 m - 1 (rij - 1)2 

The continuous eigenvalue problem requires first to consider the following symmetric case (e = 0): 

fco d^p{x, t) 



-p.p{x,t) 



("-1 



1) 



dx{^ 



(49) 



(50) 



(51) 
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with boundary condition ^a ^ — 0, p(l, t) — 0, which yields the least eigenvalue jr^^^^jw, which follows from straightforward 



algebra, see [Chapter 5] ||34(] . For the general case (mistuning design), we use the operator perturbation method [Chapter 9] ll34ll . 



the least eigenvalue for the continuous case is given by 

H=-r^^+s^ + 0{e'). (52) 

4(ni — 1)^ ni — 1 

Comparing ( |50t with ( |52] l. we have that for the symmetric case (e = 0), the eigenvalue approximation error is 0(1/ {ni — 1)^) = 

0(l/nf ), and for the mistuning design case, the error is 0(l/(ni — 1)^) + O(e^) = 0{l/n1) + O(e^). Now, take Laplace 

transform for both ( i44l i and (|45j, the characteristic equations for the coupled-ODE and PDE models are s^ + bos + A = and 

s^ + 60S + /i = respectively, which implies that the stability margin approximation error are also ©(l/nf) for symmetric 

control, and 0{l/n1) + 0{e^) for the mistuning design case. This completes the proof. ■ 

B. Simulations 

We now present results of some time-domain simulations that show the time-domain improvements - manifested in faster 
decay of initial errors - with the mistuning-based design of control gains. These simulations provide further corroboration of 
the two main conclusions of this paper: 

1) Stability margin can be improved by using a higher-dimensional information graph with symmetric control. 

2) Stability margin can be improved by using mistuned control gains for the same information graph. 

For the first set of simulations, we consider iV = 25 vehicles in a one-dimensional formation (Dg — 1). The initial position 
and velocity of each vehicle are randomly drawn from a uniform distribution on [—0.01, 0.01]. We carry out simulations for two 
distinct information graphs for the same physical formation which is consisted of 25 vehicles: a 26-node ID lattice (including 
1 reference vehicle) and 6 x 5-node 2D lattice (including 5 reference vehicles). Figure [12] (a) and (b) show the time histories 
of the relative position errors of the vehicles, for the ID and 2D information graphs, respectively. In both cases, the control 
strategy is symmetric with gains /cq = 0.01, bo — 0.05. On comparing Figure [12] (a) and (b), we see that the errors in the 
initial conditions are reduced faster with a two-dimensional information graph compared to the one-dimensional case. This 
observation is consistent with with the result of Theorem [T] 

The second set of simulations are carried out to test the effect of mistuning, for which we consider a formation with 225 
vehicles with a square 2D information graph - a 16 x 15 lattice (including 15 reference vehicles). The initial position and 
velocity of each vehicle was again chosen as a random small perturbation of the desired position and velocity. Figure [T3] (a) 
and (b) show the time history of the position errors with symmetric and mistuned control gains. For the symmetric control, 
the control gains are fco — 0.01, bo = 0.05. For the mistuning case e = 0.001, i.e., the gain fc(ij) is perturbed by ±10% 
from its nominal symmetric value fco- On comparing Figure [T3] (a) and (b), we see that the errors in the initial conditions are 
reduced faster in the mistuned case compared to the symmetric case. This improvement is consistent with what is predicted 
by Theorem [2] 

C. Disturbance propagation 

When external disturbances are present, we model the dynamics of vehicle i hy pi ~ pi — Ui + Wi, where Wi is the external 
disturbance acting on the vehicle. Each component of the disturbance is assumed to be independent. In the p, v coordinates. 
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Fig. 12. Compai'ison of symmetric control's performance in ID and 2D information graphs, (a) shows the relative position errors as a function of time for an 
25 vehicle platoon with a ID lattice (25 vehicles and 1 reference vehicles) as the information graph (see Figure [2] (a)), (b) shows the relative position en'ors 
for the same platoon (with the same initial condition) with a 2D square lattice (25 vehicles and 5 reference vehicles) as the information graph (see Figure [2] 
(b)). In both cases, the gains used are fco = 0.01 and feo = 0.05, and the initial condition is such that making the position and velocity have arbitrary and 
small (< 0.01) perturbation from the desired position and velocity. 



the closed-loop dynamics of the formation is given by 

■0 = A^ + 



(53) 



where ip :— [p^, v^]"'" is the state vector, w := [wi,W2, ■ ■ ■ ,wn]'^ is the vector of disturbances. We consider the vector of 
1^ — p, where pi ~ Pi — p*, i — 1,2, ... ,N, as the outputs: 



errors e := [pi, . 



■,PN\ 



e = Ci), C ^ [I; 0] 



The Hoc norm of the transfer function G^e from the disturbance w to the errors e is a measure of the closed-loop's sensitivity 
to external disturbance. For one-dimensional platoons, such a norm has been used previously in i7l lla,l22ll . Figure [T4l depicts 
the Hoc norm of G^e as a function of N, for the two cases described in Section IVI-BI Part (a) of the figure compares the 
Hoa norm of the one-dimensional and two-dimensional information graphs for the same formation with symmetric control. 
Part (b) of the figure compares the Hoo norm of the symmetric and mistuned control for the two-dimensional formation. The 
trends for the Hoo norm are consistent with the eigenvalue trends and the results of the time-domain simulations. In particular, 

1) The Hoa norm of Gwe is improved by using a higher-dimensional information graph with symmetric control. 

2) For a particular information graph, the Hoo norm of G^e is improved by using mistuned control gains over symmetric 
control. 

Analysis of these trends is beyond the scope of this work, and will be undertaken in future work. 



D. Other boundary conditions 

In this paper, results are derived for an arrangement of reference vehicles on one of the boundaries of a _D-dknensional 
information graph (see Assumption |2]). For a one-dimensional information graph, this means there is one reference vehicle. 
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Fig. 13. Comparison of time-domain performance between symmetric and mistuned control with the same 16 X 15 2D square information graph, (a) shows 
the relative position errors as a function of time with symmetric control; (b) shows the relative position eiTors as a function of time for the same formation 
under mistuned control with control gains chosen according to the mistuned gains with parameters /cq = 0.01, bo = 0.05, and e = 0.001 (i.e., ±10% 
variation from the symmetric gains). The initial condition is the similar to that as described in Figure 1121 Since there is a large amount of vehicles in the 
formation, for the purpose of showing the figure more clearly, only the first 10 vehicles' traces are shown. The other traces have the same trend. 



For a Z3-dimensional square information graph with N vehicles and D > 1, this means that there are No reference vehicles. 

In terms of the methodology of this paper, the arrangement of reference vehicles affects the boundary condition in the 
PDE approximation but not the PDE itself. Under Assumption |2l the boundary condition is the Dirichlet boundary condition 
at xi = 1 and Neumann boundary conditions for other boundaries of [0,1]^ (see (|28]l). More generally, the presence of 
reference vehicles on an additional boundary means that the PDE approximation will have Dirichlet boundary condition for 
these boundaries. Figure [T5] enumerates some of the possibilities for the two-dimensional case. 

It is straightforward to extend the analysis and the conclusions of the preceding sections to handle these more general 
boundary conditions. For asymptotic trend of the stability margin with N does not change with the boundary conditions. The 
presence of additional reference vehicles affects only the constant in front of the asymptotic formula. Consider for example the 
two-dimensional square information graph with N vehicles and 4\/]V — 4 reference vehicles arranged along the 4 boundaries. 
The PDE approximation is again given by JTh . The boundary conditions are all Dirichlet type. The scaling laws for this case 
are described in our conference paper iSfl. We summarize the results for the symmetric and the mistuned cases in the following 
corollaries. 

Corollary 4 (from Corollary 1 in [8]): Consider an A^-vehicle formation with dynamics ([U and control law ^ under 
Assumption [T] with nominal symmetric control gains fco and bo, whose information graph is a {^/N + 2) x {\/N + 2) 
2D lattice. Let all the nodes on each of the 4 faces of the information graph correspond to reference vehicles, so that the 
boundary conditions of the PDE ( |29] | are all of the Dirichlet type. The closed-loop stability margin is given by 



27r^J_ J_ 

^' bo N^^^N^' 



D 



Comparing this result with Corollary [T] (when D = 2) shows that the benefit of extra information (four times as many 
vehicles provided reference trajectory information) is a factor of 8 improvement in the closed-loop stability margin. 



25 



lU 










■ : 










■ 




10*^ 


■ 

• 


ID information graph 
2D information graph 


■ 
■ 




: 


10^ 




H 






, 


8 




■ 








10^ 


■ 

• 


■ 
■ 

• 

• 
• 


• 
• 


• 


• 


10^ 










- 



10 30 ,^ 100 

N 



500 



10 



o 



10' 



Symmetric 
Mistuned 



10 30 ,^ 100 

TV 



500 



(a) 



(b) 



Fig. 14. Hoo norm of the transfer function Gme from disturbance w to spacing error e. (a) compares tlie Hrx> norm as a function of N, witii ID and 
2D information grapiis, with all else remaining the same, (b) compares the Hoo norm in symmetric control and mistuned control (with ±10% mistuning), 
when the information graph is the same (a 2D lattice). In all cases, the gains used are fco = 0.01 and bo = 0.05. The mistuned gains used are those given 
in Theorem |2] with e = 0.0001. Norms are computed using the Control Systems Toolbox in MATLAB®. 
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Fig. 15. A pictorial representation of the 3 of the several possible boundary conditions for a 2D information graph. 

Corollary 5 (from Corollary 2 in [8]): Consider the same vehicle formation with the same information graph as stated in 
Corollary |4] Now consider the problem of maximizing the stability margin by designing the proportional control gains k[i^j), 
where the gains are required to satisfy \k(i^j) — fco I < e for every {i,j) G E, with e G (0, fco) being an arbitrary pre-specified 
constant. For vanishingly small values of e, the optimal control gains of the i-th vehicle (i = 1, . . . ,N) are given by the 
following formula 

fc(j ji+) = fco + 2e{H{iiCi - 0.5) - 0.5), 



«(i,ii-) 



= fco - 2e{H{iiCi - 0.5) - 0.5), 



«(i,i2+) 



fco + 2e(i7(i2C2-0.5)-0.5), 
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fc(i,i2-) = ko — 2s{H{i2C2 — 0.5) — 0.5), 

where H{x) = 1 if x > 0, H{x) = if a: < 0, and H{0) = 0.5, and ci is defined in (|22] |. The resuhing stability margin is 
given by 

o 8e 1 ^, 1 , 



The formula is asymptotic in the sense that it holds when ni,n2 —> oo and e — ?► 0. D 

Comparing this result with corollary |2] for D ^ 2 shows that with mistuning, having four times as many vehicles that have 
reference trajectory information results in a factor of 4 improvement in the stability margin. 



E. Comparison to earlier work 



9D. In 



, Bamieh et. al. proposed certain 



There are connections between the results of this paper and the results in 
macroscopic performance measures to quantify the sensitivity to disturbances of vehicular formations. The vehicles were 
modeled as double integrators and the feedback control was symmetric. The information graph considered by Bamieh et. 
al. was a D-dimensional torus, which is similar to a D-dimensional square lattice. It was shown in [9] that the measure of 



m 



disturbance amplification proposed in ||9|] grows without bound as a function of N for D = 1 and D = 2, but it is uniformly 
bounded with respect to N for D > 3. In contrast. Corollary [T] shows that there is no uniform bound on stability margin in 
any dimension for square lattices. The scaling law for the stability margin, however, improves with increasing D, as well as 
with mistuning. In summary, the asymptotic behavior of the stability margin in dimensions D = 1 and D = 2 is similar to 



that of the macroscopic performance measure of Bamieh et. al. in 
3 and higher 



. However, the trends are quite different in dimensions 



VII. Conclusion 

We studied the closed-loop stability margin with distributed control of a network of N double integrator agents. Information 
graphs (within the class of D dimensional lattices) that characterize the information exchange structure among vehicles were 
examined. We first examined the case of symmetric control, in which every vehicle uses the same control gains. For a square 
information graph, the stability margin approaches zero as 0(l/iV^/^) as A^ -> oo. Therefore, the stability margin can be 
improved by increasing the dimension of the information graph. For a non-square information graph, the stability margin can be 
made nearly independent of the number of vehicles by choosing the "aspect ratio" appropriately. The trade-off is that increasing 
the dimension of the information graph or choosing a beneficial aspect ratio may require long range communication and/or 
entail an increase in the number of reference vehicles. These results are therefore useful in investigating design trade-offs 
between performance and the cost of designing information architectures for distributed control. 

Second, a mistuning-based approach for stability margin improvement over symmetric control is proposed that consists of 
making small changes to the gains over their nominal values in the symmetric case. The scaling laws for the stability margin 
with mistuned control showed that with arbitrarily small amount of mistuning, the stability margin can be improved significantly 
over symmetric control. The mistuned control is simple to implement and therefore attractive for practical application. 

A PDF approximation was derived to aid the analysis and design that was carried out in the paper. The control design 
problem is much more tractable in the PDF domain than in the original state space domain. In particular, the PDE model 
provides insight into the effect of asymmetry in the control gains on the stability margin, which enabled the mistuning-based 
design. Such insight is difficult to gain by examination of the state-space model. Although the PDE approximation is valid 
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only for N -^ oo, numerical calculations using the PDE model show that accurate predictions are obtained even for small 
values of N. 

The information graphs studied in this paper are limited to Z?— dimensional lattices. More complex graph structures will 
be explored in future work. We believe that the PDE approximation will be beneficial here, by allowing us to sample from 
the continuous gain functions defined over a continuous domain to assign gains to spatially discrete agents. Another future 
direction of research is the examination of the closed-loop's sensitivity to external disturbances. For symmetric control, this 
issue was investigated in Igt]. Numerical tests reported in this paper show that mistuning reduces the closed-loop's sensitivity 
to external disturbances. Analysis of the effect of mistuning on the closed-loop's sensitivity to external disturbances will be 
carried out in future work. 
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Appendix 

Proof of Theorem \3] The proof proceeds by a perturbation method. Let the eigenvalues of the perturbed PDE ( |42] | and the 
Laplace transform of p{x, t) be 

respectively, where s- and 77^°^ are corresponding to the unperturbed PDE ( [3TT l. Taking a Laplace transform of both sides of 
the PDE ( |42] | with respect to t, plugging in the expressions for s and r/, and doing an 0(1) balance leads to the eigenvalue 
equation for the unperturbed PDE: 

n(0)-n wh^r^-P — l'/'J")^2 , , JO) 



Vv^"" = 0, where V := (^(s^)' + bos'f' - A 

where Co is the Laplacian operator defined in (|32] |. Recall that the solution s- ,r]'^^^ to this equation have been previously 
given. Eq. dJTj l provides the formula for si ' (i.e si), and 77^"^ = ^ (j)-j{x)aj{s), where (fi^x) is given by equation (|35T l. Next 
we do an 0(e) balance, which leads to: 

V ^^ rid - 1 Sxd ^^ 2{nd - l)^ ^x^ i i y ' 

For a solution Ty*^*^^ to exist, R must lie in the range space of the operator V. Since V is self-adjoint, its range space is orthogonal 
to its null space. Thus, we have, 

<i?,0j<f)>=O (54) 

where (t)j{x) is also the (Zi, ^2, • ■ • , IdY^ basis of the null space of operator V. We now have the following equation: 

7o "'0 ^^"rf-1 oxd ^^2{nd-iy dx^ ' z z y ' 

Following straightforward manipulations, we got: 

(60 + 2s?^)s[,") / • • • / {(bj{x))^dxi ■■■dXD = 



' ' /o Jo 



(2/1 - l)7r fi 



4(ni - 1) Jo 







k'l'lx) sin ((2^1 — l)7ra::i) cos^(Z27rx2) • • • cos^(/D7ra;_D) dxi ■ ■ ■ dxn 



hTT f^ [^ ~„,,^, _^2 / (2^1 - l)^a;i 

2(n2 - 1) io 







A;™(a;) cos" ( ) sin(2Z27rx2) • • • cos {Iottxd) dxi ■ ■ ■ dxo 



lo-n- [^ f^Jmt^x 2 ({'^h-l)'^xi 



2(nD - 1) Jo 



/ fc^(x)cos ( ) COS {l2''^X2) ■ ■ ■Six\[2li:,TTXD) dxi- ■ ■ dxD 

Jo 2 



+x •••y„ (^(;^^+---+^(^^^)'^'-<^)^^^ 

When TT-i, . . . , riD are very large, bo + 2s-. w &o- Using this, and substituting the equation above into Sj^— si +es-f +0{e^), 
we get the following: 

(0) 
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s{2h~l)n f' Tr.™,,,,, .,,, ,^_^_,2„__^ . 



Aboini - l)M Jo 







2bo{n2 - 1)M Jo 







k"^{x) sin ((2^1 — l)7ra::i) cos {l2'^X2) ■ ■ ■ cos [Idt^xd) dxidx2 ■ ■ ■ dxo 

k^{x) cos ( ) sin(2/27ra;2) cos {I^ttx^) ■ ■ ■ cos (/dttx^i) dxidx2 ■ ■ ■ dxo 



^^D-^ r r Im/^^ „„.2 ( (2^1 - 1)7^2: 



k^{x) cos ( ) • • • cos (Z(£i_i')7ra;(£)_i)) sin(2Z£)7ra;£)) dxidx2 ■ ■ ■ dxn 



2boinD~l)MJo "Jo ' 2 

e /"^ /"V ^f(^) 9277(0) kUx) 9277(0) \ ,, , ^, 

where M := J^ ■ ■ ■ J' {(j)j{x))'^dxi ■ ■ ■ dxjj = J^ ' ' ' /n '^'-'^^ ( ^~2 ) cos2(Z27i'a^2) ■ ■ • cos'^{li:i'KXB)dxi ■ ■ ■ dxjj. Without 
mistuning, the least stable eigenvalue is given by sL ^ ) with an associated eigenfunction '/>(i.o,o,...)('^) ~ cos(|-xi), which 
is almost everywhere positive in [0, 1]^. As a consequence of the Sturm-Liouville theory for the elliptic boundary value 
problems, the possibility of "eigenvalue cross-over" is precluded. That is, some other eigenvalue from becoming the least 
stable eigenvalue in the presence of mistuning is ruled out. The standard argument relies on the positivity of the eigenfunction 



corresponding to sL ^ ^ y the reader is referred to Il35ll for the details. Thus, for vanishingly small e, the least stable eigenvalue 



is 5(1.0,0,... )' even in the presence of mistuning. Setting /i = 1 and ld = Q for d > 1 in (|55] |, we obtain the result. 



